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U  •  N  I 

A  new  family  of  life  models  proposed  by  Glaser  (1980)  lias  dnrs4t  ies^ 

2 

of  the  form  C(0^  ,©2  ,0-j)cxp{  -0^t  “  ^t  +  ®3 logt],  The  family  includes  all 
gamma  distributions,  all  exponential  distributions,  all  normal  distribu¬ 
tions  left  truncated  at  zero,  and  a  variety  of  distributions  having 
bathtub  shaped  failure  rate  functions.  In  this  paper  statistical  in¬ 
ference  for  the  family  is  considered  with  an  emphasis  on  determining 
whether  the  sampled  model  fits  the  bathtub  category. 


. 


1.  INTRODUCTION 


A  family  3  of  life  distributions  has  been  proposed  by  Glaser  (1980) 
It  is  characterized  by  the  collection  of  density  functions  of  the  form 


f(t)  =  c(el,e2,e3)exp[-e1t  - ©2t2 +  e3  log  t},  o<t^<». 


(l.D 


where  C  is  a  normalizing  constant,  and  the  (natural)  parameter  space  is 
the  union  of  :  "  •  <  ^  <  ® »  ©2  >  °»  ®3  >  “  0  and 

[(03>02,03)  :  0^  >  0,  02*O,  03  >  -  1],  The  special  case  ©2  *  0  gives 
the  class  of  all  gamma  densities.  The  special  case  @2  >  0,  ©3  *  0 
gives  the  class  of  all  truncated  normal  densities  with  left  truncation 
at  0.  Usefulness  of  the  family  derives  from  its  failure  rate  proper¬ 
ties.  (Cf.  Glaser  (1980).)  I.et  r(t)  =  f(t)/F(t)  denote  the  failure 
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rate  function,  where  F(t)  =  J  f(y)dy.  If  0^  >  0  and  0^  >  0,  then  r(t) 
is  monotone  increasing  (IFR) ;  however,  if  9^  >  0  and  0^  <  0,  then  r(t) 
is  bathtub  shaped  (BT) :  precisely,  there  exists  t^  >  0  such  that  r(t) 
is  strictly  decreasing  for  t  <  tQ  and  strictly  increasing  for  t  >  tQ. 

For  the  gamma  situation  &2  =  0,  r(t)  is  monotone  decreasing  (DFR)  if 
<  0,  constant  if  0^  ■  0  (the  exponential  density  case),  and  monotone 
increasing  if  0^  >  0.  The  family  thus  constitutes  a  rich  class  of  models, 
encompassing  all  gamma  distributions  and  all  normal  distributions  left 
truncated  at  zero,  as  well  as  a  variety  of  distributions  having  bathtub 
shaped  failure  rate  functions.  The  latter  property  is  of  practical 
concern,  since  a  bathtub  failure  rate  corresponds  to  the  common  phenome¬ 
non  of  an  item's  improving  with  age  initially  and  deteriorating  with  age 
eventually.  In  this  paper  statistical  inference  for  the  family  will  be 
considered  with  an  emphasis  on  determining  whether  the  model  sampled  fits 
the  bathtub  category. 

Maximum  likelihood  estimation  of  the  parameters  based  on  random 
sampling  is  reviewed  in  Section  2.  In  Section  3  various  hypothesis 
tests  are  proposed  which  address  the  question  of  whether  a  bathtub 
classification  is  appropriate.  Finally,  in  Section  4  maximum  likelihood 
estimation  based  on  type  II  censoring  is  considered.  Included  is  an 
estimate  of  the  location  of  the  bathtub's  "plug." 

2.  MAXIMUM  LIKELIHOOD  ESTIMATION  BASED  ON  RANDOM  SAMPLING 

Maximum  likelihood  estimation  of  the  parameters  under  random  sampling 
is  discussed  by  Glaser  (1978).  The  maximization  is  complicated  by  the 


different  roles  played  by  0^  depending  on  whether  it  is  positive,  zero, 
or  negative.  A  re-parametrization  motivated  by  (ej ’s  being  basically 
a  scale  parameter  proves  effective.  The  overall  family  of  models 
comprises  the  following  four  classes: 

Class  1.  f(t) *[X/r(0,p)]exp{-Xt-0(Xt)2)(Xt)p"1,  X  >  0,0  "  0,p  '  0  (2.1) 

where  r(0,p)  is  defined  by  r(0,p)  =j*^  exp[-y  -  0y2]yP ”^dy . 

Class  2.  f(t)  *[6/A(0,p)]exp{6t -0(6t)2}(6t)p"1,  6>O,0>O,p  >0,  (2.2) 

0®  ty  | 

where  A(0,p)  is  defined  by  A(0,p) * T  expfy  -  0y  jyp -1dy . 

"0 

Class  3.  f(t) -[2X/r(p/2)]exp{-(Xt)2](Xt)p‘1,  X>0,P  >0,  (2.3) 

where  T(*)  is  the  gamma  function,  defined  by 

r(5)  ■  J*  exp[“y3y^-1dy. 

Class  4.  f(t)  =  [X/r(p)]exp[-Xt}(Xt)p_1,  X>0,p>0.  (2.4) 

Class  1  is  formed  from  0j^  >  0,  ©2  >  0;  Class  2  is  formed  from  0^0, 

@2  >  0;  Class  3  is  formed  from  0^  *  0,  02  >  0;  and  Class  4  is  formed 
from  0j  >  0,  02  *  0.  For  Classes  1,  2,  and  3,  f  is  BT  if  p  <  1,  TFR 
if  p  >  1,  and  truncated  normal  if  p  ■  l.  Class  4  is  the  class  of 
gamma  densities,  with  f  DFR  if  p  <  1,  TFR  if  p  >  1,  and  exponential 
if  p  ■  1. 

As  described  by  Glaser  (1978),  the  maximum  likelihood  estimator  of 

(0^,02 >03)  can  be  computed  by  the  following  algorithm.  Let 

t  *  (tj . tn)  denote  the  observed  random  sample.  For  each  Class  i  , 

1*1, 2, 3, 4,  let  L.(t)  denote  the  likelihood  function  based  on  the 
1  ^ 


4 


corresponding  parametr ization  (for  example,  L.(t)  *=  L.  (t  ;X  ,0,p  )  «*  n?f  ( t . ) 

1  ^  1  ^  l  j 

A 

with  f(t)  taken  from  (2.1)),  and  let  T. . (t)  denote  the  maximized  likeli- 

i  ^ 

A(l) 

hood  for  the  class.  Hence,  for  Class  1,  if  (X  ,0  ,p  )  denotes  the 

A 

MLE  of  (X,6,p)  based  on  (2.1),  then  L^(£)  is  the  maximized  likelihood  for 

Class  1,  i.e.,  L,(t)  =  L,  (t  ;X^^  ,0^\p  ^  ^)  .  The  overall  maximum  likeli- 

hood  is  L(t)  *  max,  ,  . L.(t),  and  the  MLE  of  (0, ,©„,©_)  is  the  value  of 

'v  1  <  1  <  4  1  **  LZJ 

(0  ,©2.83)  corresponding  to  the  MLE  of  the  parameters  in  the  class  with 

A  A 

the  largest  maximized  likelihood.  For  example,  if  L(t)  *  L-00 ,  then 
the  MLE  of  (©^©2,62)  is  (  -  6<2),e(2)(6(2))2,p(2)-l). 


Computational  algorithms  for  the  MLEs  for  the  individual  classes 

are  presented  by  Glaser  (1978).  To  review,  let  (s^,S2,s^)  denote 

(Je  tJf  ^  £  tj2,  ^£  log  tj).  Then  (X(1)  ,©(1)  ,p (1) )  is  the  solution 

2  2 

to  the  system  of  equations  Sj^  / 8 2  *  T  (9,p  +  l)/r(©,p  )r(0,p  +  2) , 

a 3  -  log  sL  «  t(6»P)  "  log[r(9,p  +  l)/r(0,p)],  and  X  =  r(0,p  +  l)/sir(©,p) , 

where  f(0,p)  *-8  defined  by  \|r(9,p)  *  (9/9p)  logr(0,p) 

-  [r(0,p)]“1  J*0  exp{- y -0y2}yp’1logydy.  Similarly,  (6(2) ,©(2) ,p (2)) 

2  2 

is  the  solution  to  the  system  of  equations  s^  !s^mh  (8,p+l)/A  (0,p  )A(©,p+2) , 
Sj-logSj  -  £ (0,p )-log[A (0»P+1)/A (0>P )] >  and  6 -A(0,p+l)/s1A(0,p)  ,  where 


^(0,p)  is  defined  by  £(0,p)  ■  (8/9 p)  logA(0,p) 

-  [A(e.p)]"1  J,o"exp{y-0y2}yP"1  logydy.  Also,  (X(3),p(3))  is  the 

1  l  P  1  P 

solution  to  the  system  of  equations  Sj-^logS2  ■  t  “  2  ar|d 

2 

X  ■  p/ia^,  where  f(‘)  is  the  psi  function,  defined  by 

*(?)  -  (9/9|)log  T(?)  -  Cr(F.)]-1  J*  exp(-y]yP_1  log  y  dy.  Finally,  for 

°  “(4)  “(4) 

the  well-known  gamma  parametrlzation,  (X  ,p  )  is  the  solution  to 
the  system  of  equations  s^  -  log  s^  ■  t(p)  *  log  p  and  X  ■  p/s^. 
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3.  BATHTUB  RELATED  HYPOTHESIS  TESTS 

In  this  section  assume  T^,...,Tn  constitutes  a  random  sample  taken 
from  some  distribution  in  3  .  Certain  hypothesis  tests  regarding  the 
condition  of  a  bathtub  failure  rate  function  will  be  considered. 

(a)  Testing  H  :  =  0  versus  K  :  9^  >0. 

The  subclass  of  3  restricted  to  «  0  represents  the  class  of 

gamma  distributions,  each  of  which  has  a  monotone  failure  rate  function. 

The  alternatives  02  >  0  include  bathtub  failure  rate  functions  (0^  <  0) . 

The  form  of  optimum  critical  regions  depends  on  whether  0^  is  known. 

If  0^  is  known,  a  uniformly  most  powerful  unbiased  (UMPU)  test 

2  2 

exists  based  on  the  statistic  V  *  ETj  /(ET^)  .  Since  for  the  gamma 
situation  (i.e.,  &2  -  0),  V  is  distributed  independently  of 
(cf.  Glaser  (1976)),  it  follows  from  the  general  theory  given  by  Lehmann 
(1959,  Chapter  5)  that  the  level  a  UMPU  critical  function  is 

if  v  <  C0 
if  v  >  CQ 

where  C^  is  determined  by  the  equation 

«  •  r«2-otV  <3'1' 

An  (n-l)-fold  integral  representation  of  (3.1)  involving  the  Dirichlet 
distribution  allows  computational  evaluation  of  C^,  Assume  H  :  ©2  ■  0 
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holds.  For  convenience,  denote  n-1  by  m  and  0^+1  by  p  .  Define 
Zt  -  Tj/S"  Tj»  1*1.  Then  the  var  iables  , .  .  .  .Z^  have  a  Dirichlet 


distribution  with  joint  density 


h(V...,zm>  -  [r(np)/r^(p)]n“  zP*1  (i 

where  0  <  zl  <  1,  1*1,..., m,  and  0  <S'J1  zt  <  1.  Since 


Z.2  +  (1  -  E®  Zt)2,  (3.1)  becomes 


f  •: -Jh(zi . z«n)dV--dV 


D(C0) 


(3.2) 


where  D(CQ)  -  {  (Zj, . . .  .z^)  :  0  <  z£  <  1,  1-1 . .  0  <D]L  <  1.  and 

D®  zt2+  (1  -E^z^2  <  Cq3.  Evaluation  of  CQ  by  means  of  (3.2)  is  com¬ 
plicated  by  the  awkward  limits  of  Integration.  A  tractable  form  of 
(3.2)  Is  achieved  upon  applying  an  orthogonal  transformation  which  In¬ 
corporates  Z)®  zi  as  a  variable.  Consider  the  mxm  orthogonal  matrix, 
0,  obtained  by  the  Gram-Schmidt  orthogonalizatlon  process  (cf.  Kaplan 
(1952)),  and  defined  by 


[m(m-l)]  2 


[(m-l)/m]2  -Cra(m-l)] 


■[m(ra-l)] 


-[m(m-l) 3 


0-  -0-1) (m-2)] 


C(m-2)/(m-l)]Z  -[(m-l)(m-2)]  ...  -[ (m-l)(m-2): 


■[ (o*2) (m-3)]  2 


[(m-3)/(m-2)]z  .  .  .  -[ (m-2) (m-3)] 
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The  transformation  u  *  0  z  gives 

^  j  fW  MM 

_2  "2 

u  *m  z.  +  . . .  +  m  z 
11  m 

-I  I  -I  _ I 

u2  -  -rm(m-l)]  2  Zj +  [  (m-l)/m]2  z2 -[m(m-l)]  2  z^  -  . . .  -  [m(m-l)]  2  z^ 

_  1  1  _  1  1 
u3*-[(m-l)(m-2)]  22l+[(m-2)/(«n-l)7z3  -[  (m-1)  (m-2)]  2  z^  -...-[  (m-1)  (m-2)]  2  zm 


u 


k 


-[ (m-k+2) (m-k+1) ] 


1  I  .1 

2zL  +  [ (m-k+1)/ (m-k+2) ]2zfc-[ (m-k+2) (m-k+1)]  2zk+J- 


-[ (m-k+2) (m-k+1 )]2z 


u  •  -  2 


z.  +  2  2  z  . 
1  m 


The  restrictions,  0  <  z^  <  1,  i«l,...,m,  and  0  <  1,  translate  to 

-I  -I  I 

ut  €  At,  i-1, ...  ,m,  where  Aj  -  (0,m  2)  ,  -  (-[m(m- 1)]  2,  [  (m-l)/m]2>  , 

-I  1  .11 

A3  -  (-[ (m-l) (m-2)]  2,  [ (m-2)/ (m-1)]2) . Am  -  (-  2  2,22).  The  condition 

n  2  m  2 

E”z^  +  (1  -  Z)  jZ^)  <  Cq  imposes  marginally  on  the  restriction 

I  II  1 

m2/n  -  (l/n)(nCg-  l)2  <  u^  <  m2/n  +  (l/n)(nCg-  l)2.  Consequently,  (3.2) 
assumes  the  computable  form 


(«! . um>  !ui6AinBi(Co> 


J,Cr(np)/rn(p)3Cn“z1(ul, 


)]^  ^(l-m2u.)p  ''du  ...du 
ii  1  m 


r 


(3.3) 


% 


L 


where  B^Cg)  -  (m2/n  -  (1/n) (nCQ  -  l)2,m2/n +  (1/n) (nCQ  -  1) 2) , 

1 

B^(Cq)  *  (O.Cq  -  [2  j^l^j  ^  ^  ™  u  i )  ]  ) ,  i*2  , . . . ,m,  and  z^(u^,,..,  u^) 

Is  obtained  from  z  -  O'u.  For  fixed  C„,  the  right-hand-side  of  (3.3) 
can  be  computed  by  numerical  Integration.  Since  the  right-hand-side 
increases  with  C^,  the  exact  value  of  CQ  corresponding  to  a  given  level 
a  can  be  found  by  an  easy  trial  and  error  process. 

A  large  sample  approximation  of  CQ  Is  readily  obtained  from  multi¬ 
variate  central  limit  theorems.  Assume  H  :  0j  *  0  holds  so  that 

T, , . . . ,T  constitutes  a  random  sample  from  a  gamma  distribution;  for 
1  n 

convenience,  denote  0^  + 1  by  p  and  0^  by  X  .  Since 

E0  “  <P+k-l>---P/xk  for  k-l.2.3,4, 

it  follows  from  the  multivariate  Lindeberg-Levy  central  limit  theorem 
that  under  H:0-“O,  n7(Y  -  p)  *  N,(0,E),  where 

Z  Nil  n*  Z  ^  N 

Y  isT  2  p(p  +  l)X2  2p(p  +  l)(2p +  3)/X4  2p (p  +  1) /X 

*  Q  ■  CiITj  )  •  £  ■  (  p/x  ) ■ ■  5  -  C  2p(p  +  l)/x* 

2 

By  applying  the  transformation  g(Yn)  ■  Yjn/Y2n  '  ^ol^ow8  Rao 

(1973,  p.  387))  that 
1 

n2(g(Y  )  -  g(p>>  4  N(0,  BSB’)  ,  where 

N  |1  ^ 


B  -  Og(Yn)/3Yin,3g(Yn)/aY2n)l(Y  (X2/p2,-2X(p  +  l)/p2>.  That  is, 

1 

n2(nV  -(p  +  l)/p)  £  N(0,2(p  +  l)/p3) ,  where  V  =  V  .  A  large  sample 

n  n 

approximation  of  Cq  is  therefore 


C 


0 


1  fP  +  1 
n  1  p 


where  z  is  defined  by  $(z  )  *  1  -  or .  In  computing  C_  by  the  exact 

a  a  u 

multiple  integral  representation  (3.3),  a  logical  starting  point  would 

be  this  large  sample  approximation. 

If  unknown,  the  UMPU  test  is  conditional  on  the  statistics 

2 

(T,S)  -  (ETj ,£  logTj) .  Denote  by  U.  From  the  general  theory  of 

Lehmann  (1959  ,  Chapter  4),  the  critical  function  of  the  level  a  UMPU 
test  is 


Cp(u,t,s 


1 

0 


if  u  <  CQ(t,s) 

if  u  >  c0(t»s) 


where  CQ(t,s)  is  determined  by  the  equation 


a  -  P^^U  <  C0(t,s)|(T,S)  -  (t,s)}.  (3.4} 

Consider  the  conditional  distribution  of  (T. . T  ),  given  (T,S)  *  (t,s), 

i  n 

assuming  H  :  0^  ■  0  holds.  Again  for  convenience  denote  0^  by  X  and 
0^  +  1  by  p  .  The  joint  density  of  (Tj,...,Tn)  is  [Xnp/rn(p)]e^p"1^8"^t, 


in 


Because  of  the  independence  of  and  (e^n/j^T)n  ,  it  follows  from 
Glaser  (1976)  that  the  joint  density  of  (T,S)  is 

[Xnp/r"(p)][(2n)<"-1>/2n1/2/r(2ji)]t-V8-Xt(n  logi.S)<n-3>/2fn(o’/(V), 

for  a  certain  function  §n(*).  Consequently,  the  conditional  density  of 
<T1,...,Tn),  given  (T.S)  =  (t.s),  is 

r_,n-l.  ...  .  (n-l)/2  1/2,..  -s,  ,  t  .-(n-3)/2  -1.  s..t,n. 

[r(— )/(2n)v  n  jte  (niog--s)  (<?/(“)  ).  The 

essential  point  is  that  the  conditional  density  is  constant  (i.e., 

uniform)  in  (tj,...,tn)  over  the  region  for  which  £tj  »  t  and  Ilogt^*s. 

The  right-hand-side  of  (3.4)  is  thus  the  ratio  of  (n-2)-dimens ional 

volumes , 


Volume{  (tlf. . .  ,tn)  :  Et^2  <CQ(t,s)  .Et^-t.E  log  t^-s] 
Volume{(tlt . . . ,tn>  :  Et^-t.E  log  t^-s) 


Unfortunately,  computation  of  these  volumes  appears  prohibitively  com¬ 
plicated.  An  exception  is  the  rather  uninteresting  case  n=3,  where  a 
cubic  equation  can  be  obtained  which  generates  C^. 

A  large  sample  approximation  for  Cg(t,s)  can  be  found  by  methods 
similar  to  those  used  in  the  known  0^  case.  Let 

Y  «  (—  ET  2  —  ET.,  —  ElogT.)'.  It  follows  from  the  multivariate 

~n  vn  j  n  j’  n  y 

Lindeberg-Levy  central  limit  theorem  that  under  H  :  9^  m  0, 

1 

n2(Y  -u)  £  N_(0  ,  E  ) , 

Mil  M  J  M  W 


I-  « 


p(p+l)/x' 

where  p,  p/X 

iJ((p)-log  X' 


and  H  = 


2p  (p+1)  (2p+3)/X 
2p(p+l)/X3 
(2p+l)/X2 


2P(P  +  1)/X 
p/X2 
1/X 


(2p  +  l)/>.  , 
1/X 

v'OO  ' 


More  importantly,  the  relevant  conditional  distribution  converges  to 

normality.  Partition  Y  ,  p,  and  T.  as 

~n  ^  ^ 


Y  - 

/-u  \ 

n  n 

»  H  * 

f"  } 

X  p(p+i)/\2  X 

/Ml 

-T 
n  n 

V-S  > 
x  n  n 

rm 

,  (2) 
vJi  ✓ 

p/X 

^♦(pj-log  X  ) 

rhi  i  EuA 

1 

X2p(p+l)(2p+3)/X4 

2P(p+l)X3 

(2p+l ) /X2\ 

I* 

**■* 

1 

f 

1 

1 

1 

• 

- 

2p(p+l)/X3 

p/X2 

1/X 

^  ~2i  ;  £22  j 

^  (Zp+1)/X2 

1/X 

*'(p)  y 

From  the  work  of  Steck  (1959),  it  follows  that  the  limiting  conditional 


l  -T 
n  n. 


distribution  of  n2(^Un  -  p^),  given  n2^  ^  ^  -  j^2^ 


-S 
n  n 


1  ~  t 

"Z  ss-n  n 


“  n2(Cl  )  •  £(2))*  l*e”  glven  (sn)  *  Cs")’  18  univariate  normal 


-  s 
n  n 


,  - 1 
-1  ,,n  n 


2  //  \  (2)^  2  -1 
with  mean  H  "  "  £12  S22  CC  1  )’£  )  a"d  varlance  ®  -Sn“Sl2  S22  hi  ' 


-  s 
n  n 


5 

J 


m 


!SSm 
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The  difficulty  introduced  by  the  fact  that  the  asymptotic  conditional 
distribution  depends  on  the  unknown  parameters  p  and  X  will  be  overcome 
by  substituting  estimators  which  are  consistent  under  H  .  As  shown  by 
Neyman  (1955),  such  substitution  does  not  affect  the  asymptotic  proper¬ 
ties  of  the  test.  When  H  :  ©2  *  0  is  true,  the  (consistent)  MLE,  (p,X), 


-  t 
n  n 


of  (p,X)  is  the  simultaneous  solution  to  the  equation  (  'j  =  p^ 

A  o  '  ~ 


—  s 
n  n 


(cf.  Class  4  of  Section  2).  Consequently,  the  (consistent)  MLE,  p,  of 

A  2 

p  under  H  is  p  *  0.  On  the  other  hand,  computation  of  a  yields 

a2  ■  (p/X^)[2(p+l)  -  (p(/ '  (p)-l)"1],  so  that  the  (consistent)  MLE,  a2  , 

of  a2  under  H  is  7=  (p /X4)  T  2  (p+1)  -  (p*  '  (p  )  -  l)"1]  .  Finally,  the 

(consistent)  MLE,  p^,  of  p^  under  H  is  ^p^  *=  p(p+l)/X2.  To 

summarize,  if  H  :  9^  *  0  is  true,  the  approximate  conditional  distribu- 


I  /\  A  T  t 

tion  of  n2(„un  "  )/a,  given  is  norn,a^  with  mean  0 

n  n 

and  variance  1.  A  large  sample  approximation  of  the  level  or  critical 

7 

value  Cq  ( t ,  s )  ,  from  (3.4),  is  therefore  CQ(t,s)  ■  hp  -n  z  o  . 


(b)  Testing  II  :  >  0  versus  K  :  9j  <  0  . 

If  02  >  0,  the  subclass  of  J  restricted  to  0^  >  0  has  only  IFR 
fuhetions,  whereas  the  subclass  categorized  by  0^  <0  has  exclusively 
BT  functions.  The  test  of  H  versuB  K,  assuming  0£  >  0,  is  therefore 
a  test  of  increasing  failure  rate  against  bathtub  alternatives.  For 
the  case  @2  *  0,  i.e,,  the  gamma  distribution,  the  test  of  H  versus  K 
is  equivalent  to  testing  IFR  versus  DFR.  A  UMPU  test  for  this  situation 


n 


is  provided  by  Glaser  (1973).  For  the  remainder  of  this  section  it  will 

be  assumed  that  0^  is  positive  but  unknown. 

2 

Based  on  the  notation  (U,T,S)  *  (E T ^  ,E  ,£  log T  ) ,  the  level  cr 
UMPU  test  is  conditional  on  (T,U): 

if  s  <  CQ(t,u) 
if  s  >  CQ(t,u) 

where  Cg(t,u)  is  determined  by  the  equation 

or  -  pe3=o^s  (3.3) 


cp(s,t,u)  -  | 


When  ©3  »  0,  the  sample  T^,...,Tn  is  truncated  normal;  however,  the  con¬ 
ditional  distribution  of  S  given  (T,U)  *  (t,u)  does  not  appear  tractable. 
Fortunately,  a  large  sample  approximation  of  Cg(t,u)  can  be  found  by 
employing  the  method  used  in  testing  ®  with  0^  unknown,  hot 


n  n 

-£T. 

as 

-  T 

n  j 

n  n 

v  i  U  J 

V  n  n  ' 

From  the  multivariate  Lindeberg- 


Levy  central  limit  theorem  it  follows  that  if  0^  ■  0,  then 

1 

2  £ 

n  (Y  -  u)  _  N,(0  ,  £),  where  u  is  the  mean  vector  and  £  the  covariance 

2 

matrix  of  (log  T  ,T  ,T  )',  assuming  03»O.  Partition  Y  as 
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Y 


n  Sn 


-  T 
n  n 


x  n  n 


and  partition  u  and  £  accordingly.  Then  (cf.  Steck 


(1959)),  for  0  =  0  the  relevant  asymptotic  conditional  distribution  is 
3  l  . 

normal:  i.e.,  the  conditional  distribution  of  n^(— S  -  u^)/o,  given 

T  tn 

(  )  -  (u  ),  tends  to  N(0,1),  where  ^  and  c  denote  the  respective 

n  n 

(consistent)  MLEs  of  and  a2  *  “~12  £22  ^21*  comPutet*  under  the 

assumption  that  0_  ■  0.  The  large  sample  approximation  of  Cn(t,u),  from 
J  l  u 

ill  7  * 

(3.5),  is  therefore  Cg(t,n)  =  n  p,  '  -  nz  a  . 

xfx  *  ^  ^ 

Computation  of  pv  '  and  a  is  somewhat  tedious.  The  MLE,  (0^,02), 

of  (0^,02)  is  the  simultaneous  solution  to  the  equation 


(2)  _ 


EVoTj 

Ee3-oTj 


(  C  -  01)/2  @2 
[l-0.(C-01)/20,]/20, 


where,  from  (1.1),  C  -  C^.O^O)  *  {(n/02)2exp(012/402)[l  -  «  (0t(202)  2)]) 

Suppose  for  definiteness  that  (  )  leads  to  a  positive  value  for  0^.  By 

n 

the  invariance  property  of  maximum  likelihood,  the  MLEs,  p.  and  a  ,  are 

obtainable  from  the  parametrize Cion  (inspired  by  Class  1  of  Section  2), 

2  2 
[X/r(e,l)]exp{-Xt -0(Xt)  ],  X  >  0,  ©  >  O.  Here  X  -  eL  and  8  -  0^  , 

*  »  *  *  2 

so  that  X  -  0j  and  0  ■  02^ei  *  since  this  parametrization  gives 
p*1*  ■  JJ’[X/r(e,l)]exp(-Xt -0(Xt)23logtdt  -  <i(0,l) -log  X,  and 
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f 


t’«M) 


_L  r.L9,3). 


x  nia>t*<9’2,‘*<9’l)] 


V 


;2 


_L  r(e.i)r(9.3)-r  (9.2) 

2  2 

x  r  (9,1) 

I  r(9.i)r(9.4)-r(9.2)r(9.3) 

X3  1^(9,!) 


1  r(9,i)r(».4)-r(9.2)r(9.3) 
x3  r2(9,i) 

I  r(9>i)r(9,5)-r2(9,3) 

4  2 

x  r  (9,i) 


(1)  *  * 

it  follows  that  ^  '  ■  f(9,l)  -  log  \  and,  after  a  bit  of  computation, 

*  *  2  2  2 

a  -  Ct’(9,l)  -  [g2  (p^p^  (S^j-gj  )  -  Zg^^-pp  (pj-p^  (g^-g^) 

2,  .2,  2. n  /r  2..  2.  .  .2.,  1/2 

+  83  (P3~Pj)  ^81^83*82  Jj/LS^-gj  )(8185-83  )  ~  (gl84-g2g3)  ]}  ,  where 

A  A  * 

gt  ■  r(9,i)  and  p1  (9,1).  Analogous  computations  of  p  and  a  are 

A  *  A 

possible  based  on  0^  negative  or  0^  zero.  For  negative  0^  a  parametriza- 

A 

tlon  of  the  Class  2  variety  would  be  used;  for  0^  -  0,  a  Class  3  parametri¬ 
za  t  ion  would  be  used. 


4 .  CENSORED  SAMPLING 

(a)  Maximum  likelihood  estimation  of  the  parameters  based  on 
type  II  censored  sampling. 

A  common  sampling  technique  in  life  testing  is  type  II  censoring  at 

r  out  of  n:  n  items  are  simultaneously  put  on  test;  the  lifetimes  of 

only  the  first  r  items  to  fail  are  recorded.  In  effect,  only  the  first 

r  order  statistics,  <  ...  <  T^,  are  observed  from  the  random  sample 

T,,...,T  ,  of  n  lifetimes.  To  estimate  the  parameters,  (9, ,0„,9,),  of 
in  l  L  i 

the  sampled  distribution  in  3  by  the  method  of  maximum  likelihood,  the 
basic  strategy  described  in  Section  2  for  random  sampling  can  be  used; 
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l.e.,  likelihood  maximization  is  carried  out  individually  £or  each  of  the 
four  classes,  whereupon  the  largest  of  the  four  maxima  provides  the  MLE 

AAA 

(e^,©2>®3)*  For  t*ie  *a'te  brevity,  only  Class  1  will  be  considered. 
Maximization  for  the  other  classes  is  analogous,  and  in  the  case  of 
Class  4  (the  gamma  distribution),  already  worked  out  by  Wilk  et  al. 
(1962). 

For  convenience,  denote  by  Yj  the  Jch  order  statistic  T^.  The 
(Class  1)  likelihood  function  is 


n! 


(n-r)! 


rr(e,P) 


exp 


{.fxv.SaV’}[na,J)]’  { j£ ?J^*p{-u-e<u)2}at>‘,'ldt]' 


Following  the  general  approach  of  Wilk  et  al.  (1962),  define  5  »  \yf  , 

it  r  x 

P  -  (II  yi)1/r/y_  ,  8  -  S  y./ry  ,  and  t  -S  y.2/ry  2  .  Then 

1  J  J  *  1  J  * 


log  L  "  log!"  — — — 

L  (n-r)  I  y/ 

2 

+  np  log  5  -  n  logr(©,p)  +  (rp-r)  log  p-rfs  -  r  9  f  t 
+  (n-r)  log  J (W,p,5>. 

where  J(S,p,?)  ■  J  x^  exp{  -  $x-9(5x)  )dx  .  The  equations  (d/dp)  logL-0, 

(d/d®)  logL  -  0,  and  (d/dt)  logL  •  0  are  equivalent,  respectively,  to 
the  equations 
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p"exp{ine7p)  ^  r(9,p)  los  5-  <n'r>(j^  J(«.P.?))/J(e,p.t)])  (*.D 

t*(r52)  [neTpT  a©  r(©,p)  +  (n-r)F;2j(e,p  +  2,F;)/j(e,p,«-.)]  (4.2> 

■t’1[  2r0|t-np/|  +  (n-r)  [J(0,p  +  1 ,  %)  +  2  0  5  J  (0,p  +  2.  ?.)]/J(0,p  ,  ?)]  ,  (4.3) 


where  (3/dp)r(0,p)  -  J*  exp{-y  -  0y25yp-1  log  y  dy  and  (3/3p)J(0,p , F) 

0 

00  J  2 

=  J  xp  exp{- 5  x  -  0(§  x)  }  log  x  dx.  The  MLEs  of  p,0,  and  £  are  defined 
by  the  simultaneous  solution  of  the  equations  (4.1),  (4.2)  and  (4.3). 
(Note  that  the  MLE  of  \  is  the  MLE  of  5  divided  by  y^.)  However,  an 
inverse  iteration  scheme  for  determining  the  MLEs  is  more  convenient: 
for  any  assigned  p,  0,  and  5,  use  (4.1)  to  calculate  p,  (4.2)  to  calcu¬ 
late  t  ,  and  (4.3)  with  t  to  calculate  s  .  Continue  with  the  Iteration 
process  until  values  of  (p ,  0,  §)  are  assigned  which  give  resultant 
(p,t,s)  calculations  equal  (to  a  desired  degree  of  accuracy)  to  the 
actual  observed  values. 


(b)  Maximum  likelihood  estimation  of  the  plug. 

For  a  model  in  3  having  a  bathtub  shaped  failure  rate  function  a 
parameter  of  Interest  is  the  point,  t,  known  as  the  plug,  where  the  failure 
rate  function  takes  on  its  minimum  value,  i.e.,  the  point  of  transition 
which  separates  decreasing  failure  rate  from  increasing  failure  rate.  Let 
g( • )  denote  the  reciprocal  of  the  failure  rate  function.  Then,  as  shown 
by  Glaser  (1980), 


If  ©2  >  0  and  ©3  <  0.  the  model  has  a  bathtub  shaped  failure  rate  function, 
and  the  plug  t  is  the  unique  value  of  t  satisfying  g'(t)  *  0.  Since  in 
this  case,  g'(t)  >  0  for  t  <  t  and  g'(t)  <  0  for  t  >  t,  computation  of 


T  from  (4,4)  may  proceed  by  a  swift  trial  and  error  scheme.  Similarly, 

A  A 

if  the  MLEs ,  02  and  0^,  obtained  either  by  random  sampling  or  censored 

A  A  A 

sampling,  satisfy  02  >  0  and  @3  <  0,  then  the  MLE,  T ,  of  T  is  the 

AAA 

solution  to  g' (t)  “  0,  where  (4.4)  is  used  with  (0^,02>03)  substituted 
for  (^VV* 


r 
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